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We study quenching dynamics of a one-dimensional transverse Ising chain with nearest neighbor 
antiferromagentic interactions in the presence of a longitudinal field which renders the model non- 
integrable. The dynamics of the spin chain is studied following a slow (characterized by a rate) 
or sudden quenches of the longitudinal field; the residual energy, as obtained numerically using 
a t-DMRG scheme, is found to satisfy analytically predicted scaling relations in both the cases. 
However, analyzing the temporal evolution of the Loschmidt overlap, we find different possibilities of 
the presence (or absence) of dynamical phase transitions (DPTs) manifested in the non-analyticities 
of the rate function. Even though the model is non-integrable, there are periodic occurrences of 
DPTs when the system is slowly ramped across the quantum critical point (QCP) as opposed to 
the ferromagnetic (FM) version of the model; this numerical finding is qualitatively explained by 
mapping the original model to an effective integrable spin model which is appropriate for describing 
such slow quenches. Furthermore, concerning the sudden quenches, our numerical results show that 
in some cases, DPTs can be present even when the spin chain is quenched within the same phase 
or even to the QCP while in some other situations they completely disappear even after quenching 
across the QCP. These observations lead us to the conclusion that it is the change in the nature of 
the ground state that determines the presence of DPTs following a sudden quench. 

PACS numbers: 75.10.Jm, 05.70.Jk, 64.60.Ht 


Following the remarkable advancement of the experi¬ 
mental studies of ultracold atoms trapped in optical lat¬ 
tices mm, there is a recent upsurge in the studies of non¬ 
equilibrium dynamics of closed quantum systems, in par¬ 
ticular from the viewpoint of quantum quenches across a 
quantum critical point (QCP)[3l|4j. The relaxation time 
of the quantum system diverges at the QCP resulting in 
a non-adiabatic dynamics and proliferation of topological 
defects in the final state reached after the quench. 

According to the Kibble-Zurek (KZ) scaling relation 
[3E]j generalized to quantum critical systems mm, when 
a d-dimensional quantum system, initially prepared in 
its ground state, is driven across an isolated QCP, by 
changing a parameter of the Hamiltonian in a linear fash¬ 
ion as t/r, the density of defect satisfies the KZ scaling 
^-du/(zu-\-i).^ here, v and z are the correlation length and 
the dynamical exponent associated with the QCP respec¬ 
tively [aHii]. Subsequently several modifications of the 
scaling have been proposed mMM- Similarly when the 
system is quenched to the gapless QCP, the residual en¬ 
ergy (the excess energy over the ground state of the final 
Hamiltonian) scales as on the contrary, 

when quenched to the gapped phase, the residual en¬ 
ergy follows a scaling relation identical to that of the 
defect density. Similar scaling relations for the residual 
energy and the defect density have also been derived us¬ 
ing an adiabatic perturbation theory for a sudden quench 
of small magnitude [15]. (See review articles [TMl8] ). 

It is well established that the phase transition in a ther¬ 
modynamic system is marked by the non-analyticities 
in the free-energy density whose information can be ob¬ 
tained by analyzing the zeros of the partition function 


in a complex temperature plane as proposed by Fisher 
m- These zeros of the partition function coalesce into a 
line (or area [5D]) in complex temperature plane, crossing 
the real axis in the thermodynamic limit; these crossings 
mark the non-analyticities in the free-energy density. A 
similar observation was made earlier by Lee-Yang [21j for 
a complex magnetic plane. In a similar spirit, a recent 
work by Heyl et al. |22j introduced the notion of dynam¬ 
ical phase transitions (DPTs) in connection to quantum 
quenches probing the non-analyticities in the dynamical 
free energy in the complex time plane. The idea stems 
from the similarity between the canonical partition func¬ 
tion 

Z{P) = Tr e-^^, (1) 

of an equilibrium system (where /3 is the inverse temper¬ 
ature) and that of the overlap amplitude (the Loschmidt 
overlap (LO)) defined at an instant of time t as 

G{t) = (^o|e-*^‘|V'o), (2) 

where, in the above equation, H is the final Hamilto¬ 
nian of the system reached through a sudden quenching 
of parameters, while j'i/'o) is the ground state of the ini¬ 
tial Hamiltonian. Generalizing to the complex time (z) 
plane, one can define the dynamical free energy, /(z) = 
— In G(z); one then looks for the zeros of the G(z), known 
as Fisher zeros, and can claim the occurrence of DPTs (at 
real times) when the lines of Fisher zeros cross the imag¬ 
inary axis. These DPTs are manifested in sharp non- 
analyticities in the rate function (/(t) = — In |G(t)p/A^) 
at those instants of time. This usually happens when the 
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system is quenched across the QCP |23]. The initial ob¬ 
servation by Heyl et al. [22] for a transverse Ising chain 
led to a series of works for both integrable and the non- 
integrable spin chains where DPTs were observed 

for sudden quenches across the QCP. Although, a later 
work m showed that DPTs can occur even when the 
system is quenched within the same phase. These stud¬ 
ies have also been generalized to two-dimensions |28l |29] 
where topology may play a non-trivial role |28j . 

A pertinent question at this point is how does the DPT 
depend on the integrability of the model under considera¬ 
tion or the nature of driving (slow or sudden)? Is quench¬ 
ing across a QCP essential to observe this? In this let¬ 
ter, we shall address these issues in the context of a spe¬ 
cific non-integrable model. We note in passing that the 
quantity |G(t)p denotes the Loschmidt echo which has 
been studied in recent years in the context of decoherence 
both in equilibrium [30II33j and non-equilibrium situa¬ 
tions [34U37] and has been generalized to finite tempera¬ 
ture |3H| • The LO is also connected to the work-statistics 
[39| and the entropy generation following a quench |40j . 

The model we consider here is a one-dimensional Ising 
model with a nearest neighbor antiferromagnetic (AFM) 
interaction J (scaled to unity in the subsequent discus¬ 
sion) subjected to a transverse field (P) as well as a lon¬ 
gitudinal field (h). It is described by the Hamiltonian 

m 

H = Y^ uf af+i - P ^ af - ft ^ af, (3) 

i i i 

where aiS are Pauli matrices. Por ft = 0, the model is in¬ 
tegrable with QCPs at P = Pc = ±1, while any non-zero 
value of ft renders the model non-integrable. Purther- 
more, since the APM interaction and the field ft compete 
with each other, there is a QPT from APM ordered phase 
to the disordered phase at a particular value of Pc(ft) for 
a given value of ft. As a result, one finds a phase di¬ 
agram in the P — ft plane (separating the ordered from 
the disordered paramagnetic phase, see the supplemen¬ 
tary material (SM)) starting from the integrable QCP 
(at Pc = 1, ft = 0) at one end and terminating at first 
order transition points at P = 0, ft = ±2 on the ft-axis. 

In this letter, we shall restrict our attention to the 
case when P is fixed to Pc = 1 so that the system is at 
the QCP when ft = 0; ft is driven slowly (i.e., defined 
by a rate r“^) or suddenly in the vicinity of the QCP 
when the system is always initially prepared in its ground 
state. In the presence of a small ft, a gap (AA) opens 
up in the energy spectrum and a perturbation theoretic 
calculation, valid for small ft, yields ^ ft^'* 
where Vh is the correlation length exponent associated 
with the relevant perturbation ft and z is the dynamical 
exponent associated with the QCP at P = 1. Noting 
that z = 1, one concludes that the exponent = 2. We 
note that a similar study was reported in reference |42] 
for a ferromagnetic (FM) Ising chain in a skewed field 
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FIG. 1: (Color online) Scaling of residual energy for slow and 
sudden quenching by changing the longitudinal field h with 
the transverse field F = Fc = 1, as obtained from DMRG 
studies. When ft is changed linearly as —t/r to the QCP 
(ft = 0), Eres ~ which is in perfect agreement with the 

KZ Scaling. In the Inset, we verify the scaling of Eres ~ ft^ for 
a sudden quench starting from the QCP. 



FIG. 2: (Color online) Our numerical results show prominent 
periodic occurrences of DPTs when the longitudinal field is 
slowly ramped from a large positive value to a large negative 
value as ft ~ —f/r. This periodic pattern can be qualita¬ 
tively explained by studying 7(f) of the equivalent integrable 
Hamiltonian Q as demonstrated in the SM. 

(having both P and ft); however, the universal behavior 
associated with the FM case is different. 

Our results establish that for a sudden quench start¬ 
ing from the QCP as well as a slow quench up to the 
QCP, numerically obtained residual energies per spin ex¬ 
hibit scaling relations which perfectly match earlier pre¬ 
dictions. On the contrary, there is a series of interest¬ 
ing and unexpected results concerning the scenario of 
DPTs following these quenches which are not reported 
before. Even though the model is non-integrable, we 
find prominent existence of DPTs when the longitudi¬ 
nal field is slowly ramped across the QCP. This is re¬ 
markable, given the fact that in the FM case [H] sharp 
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FIG. 3: (Color online) Numerically obtained l{t) showing 
the absence and occurrence of the non-analyticities (DPTs) 
in different situations when h is suddenly quenched, (a) No 
DPTs are observed for a sudden quench of small amplitude of 
h even if the system crosses the QCP in the process; (b) DPTs 
occur when h is quenched within the same phase; (c) DPTs 
also appear when h is quenched from a large positive value to 
the QCP at h = 0; (d) A regular (but not periodic) occurrence 
of DPTs is observed when h is suddenly quenched from a large 
positive to a large negative value. The inset shows that the 
DPTs are rounded off when the quench amplitude is even 
larger, leading to Rabi osscillations. 


the adiabatic perturbation theory |15j , for such a sudden 
quench of small magnitude starting from the QCP, Cres 
should scale as h'^h{d+z)^ gg exponent does not 

exceed 2; this is indeed true in the present case and as a 
result the exponent saturates to 2. 

Having established the scaling of Cres for both slow and 
sudden quenches, we now probe the scenario of possible 
DPTs. The Loschmidt overlap at an instant t (where the 
initial time t = 0 is set immediately after the quench¬ 
ing is complete) is given by G[t) = (^/>/| exp(—ii7^t)|'!/)j) 
and consequently one defines the rate function I{t) = 
— In |G'(t)p/7V, and investigates its temporal evolution 
to probe the signature of possible DPTs (namely, the 
non-analyticities in I(t))). Results obtained for the slow 
and sudden quenches obtained by using t-DMRG are pre¬ 
sented in Figs. I^an we below analyze the remarkable 
findings. 

We first analyze the slow-quenching of the model ([^ 
with h ~ —t/r, where h is varied from a large positive 
to a large negative value. Referring to the Fig. we 
find that I{t) shows non-analyticities which appear at 
regular (and periodic) intervals when r 3> 1 in contrast 
to the FM case [12]. To analyze this, we recall that for a 
sufficiently slow driving the dynamics is always adiabatic 
except in the vicinity of the QCP {h <C 1) where the 
relaxation time diverges. Remarkably, the non-integrable 
Hamiltonian ([^ can be mapped to an effective integrable 
model for h 1, described the Hamiltonian: 


non-analyticities are present in / (t) in the integrable case 
for /i = 0 when F is ramped across the QCP; on the con¬ 
trary, those get smoothened out when the skewed field 
is quenched through the QCP at P = 1 so that the sys¬ 
tem is always non-integrable except at the QCP. On the 
other hand, for sudden quenches the DPTs are found to 
occur whenever there is a difference in the nature of the 
ground states of the initial and the final Hamiltonians 
irrespective of the fact whether the system is quenched 
across a QCP or not. 

Let us first consider the situation when the field h is 
ramped linearly to the QCP {h = 0) as h = —ijr fixing 
P = 1. Denoting the final Hamiltonian Hf with ground 
state energy E^, and final wave function of the system (of 
length N) reached after the quench as |^/’/), the residual 

energy per spin is defined by Eres = {{'4’f\Hf\tpf)—E^)/N. 

Using the t-DMRG calculations with an open boundary 
condition, we find Eres ~ (see main part of the 

Fig. 0 ). This is in perfect agreement with the KZ scal¬ 
ing prediction, Eres ~ T-’^id+z)/Gz+i) = 2 

and z = 1 (see the discussion in the SM). We now turn 
our attention to the sudden quench, in which the system 
is initially at the QCP and suddenly a small longitudinal 
field h is switched on; in this case, numerically we find 
Eres ~ (Inset, Fig. 0). According to the prediction of 


TLeff = (1 - bh^) ^ ^ rf, (4) 

i i 

where r^’s are Pauli spin matrices and 6 is a constant 
which in our case can be chosen to be of the order of 
unity; this mapping to the model 0 is shown to ex¬ 
actly describe the low-lying excitations of the Hamilto¬ 
nian 0 in the thermodynamic limit |41j . Consequently 
so far as the slow quenching is concerned, when the dy¬ 
namics is non-adiabatic only in the vicinity of a QCP, 
one can work with the effective Hamiltonian 0 which 
represents an AFM transverse Ising chain and is equiv¬ 
alent to a FM transverse Ising chain by a simple gauge 
transformation. Both the models are exactly solvable 
by Jordan-Wigner transformation. Focussing only at the 
QCP at h = 0 and considering a slow ramp of h from 
a large positive value to a large negative value with the 
system initially in its ground state, one can derive the 
final wave function by numerically integrating the cor¬ 
responding Schrodinger equation; the rate function thus 
obtained indeed shows occurrences of the DPTs thereby 
qualitatively explaining the phenomena we observe here 
(see the SM for details). 

Interestingly, the mapping to the effective Hamiltonian 
0 also enables us to explain the absence of DPTs follow¬ 
ing a sudden quench of small amplitude across the QCP 
of the original model as presented above in Fig. 0a) be- 
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cause the interaction term in the equivalent Hamiltonian 
Q does not change sign which implies that this quench¬ 
ing does not take the system across the QCP of Hamilto¬ 
nian Q- This explains the absence of DPT in this case 
though there is a crossing of QCP in the original model. 
Though the mapping to the equivalent Hamiltonian is 
strictly valid for /i 1, in Fig. [^a), we show this ar¬ 
gument can be extended to explain the absence of DPTs 
when h is quenched from +0.7 to —0.7 crossing the QCP 
at h = 0. 

Analyzing the original Hamiltonian (§ , we note that 
the ground state is paramagnetic with all spins polarized 
in the direction of h, when /i 3> 1; on the contrary, it is a 
quantum paramagnet with majority of spins orienting in 
the direction of P when /i <C 1. The change in the nature 
of the ground state is reflected in DPT, irrespective of the 
fact whether the system crosses the QCP in the process 
of quenching. In Fig. j^b), we find a prominent presence 
of DPTs when h is quenched from 3 to 0.2; here, even 
though the quenching does not take the original Hamil¬ 
tonian across a QCP, the nature of ground state changes. 
Similar DPTs are observed when quenched to the QCP 
also (Fig. I^c)). No such DPT is found to occur when 
the nature of the ground state is the same (e.g., when 
h is changed from 3 to 2). Finally, when h is suddenly 
quenched from -|-3 to —3 across the QCP, one finds a reg¬ 
ular (but not periodic as shown in Fig. occurrence of 
DPTs (see Fig.j^c)). This is a generic feature of a sudden 
quench across the QCP as also observed in the FM case 
[24] (while the periodic pattern is only a characteristic of 
the integrability of the underlying Hamiltonian). In this 
case, the initial and final ground states are nearly fully 
polarized states with their overlap being exponentially 
small with the system size; this difference of the ground 
states results in observed DPTs. When the quench am¬ 
plitude is further increased (e.g., h = -|-5 to —5; see the 
inset of Fig. [^d)), both the initial and final Hamiltoni¬ 
ans essentially reduce to an assembly of non-interacting 
spins; in such situations DPTs are rounded off leading to 
Rabi oscillations between two fully polarized states. 

Finally, we summarize the results: we have estab¬ 
lished that Eres satisfies universal scaling relations for 
both sudden and slow quenches. Furthermore, for the 
slow quenches, the model (§ provides a unique exam¬ 
ple where one can work with an equivalent integrable 
model for r ^ 1. This mapping enables us to explain the 
KZ scaling and also a periodic occurrence of DPTs for a 
slow quenching across the QCP. This is remarkable in the 
sense that, to the best of our knowledge, the presence of 
DPTs following a slow quench of a non-integrable model 
has not been reported earlier; in the FM situation, these 
non-analyticities get smoothened out [42]. Concerning 
the sudden quench, we also present some remarkable ob¬ 
servations: in some cases, DPTs do not occur even when 
the system is quenched across the QCP; but they may ap¬ 
pear when the system is quenched within the same phase 


(even to the QCP). For very large amplitude quench of 
h across h = 0, DPTs get rounded off. These observa¬ 
tions lead us to the conclusion that concerning the sud¬ 
den quenches, it is the change in the nature of the ground 
state that is responsible for DPTs. 

We would like to conclude with the note that the 
Hamiltonian ^ has been experimentally studied using 
Bose atoms in an optical-lattice |23|, with F <C h. The 
field F of the equivalent spin chain is determined by the 
hopping amplitude t of the Bose atoms and is given by 
23/2^; F is necessarily kept small to stablize the Mott 
state necessary for the realization of a spin system. On 
the other hand, a quantum Monte-Carlo study jH] shows 
that in one dimension it should be possible to achieve a 
field F « 1. Therefore it should be possible to verify 
some of the situations of the present study in experimen¬ 
tal systems. 
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In this supplementary material, we shall illustrate how to derive the Kibble-Zurek scaling relation of the residual energy of 
the original Hamiltonian 

H = '^a4a4+i (SI) 

i i i 

when the longitudinal field h is quenched as —ijr, and also show how the behavior of the Fisher zeros in the complex time 
plane dictates the occurrence of DPTs for such a slow quench. The transverse field F is always set to unity so that the system 
is at the integrable quantum critical point (QCP) when h = 0. The phase diagram of the model and the quenching path are 
shown in Fig. |S1| 


h 
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FIG. SI: The schematic phase diagram of the model given in Eq. (Sll; the solid line extending from (F = 0, h = 2) to 
(F — 0,h = —2) through (F = 1, /i = 0) separates the antiferromagnetic (AFM) phase from the paramagnetic (PM) phase. The 
points (F = 0, h = ±2) denotes the first order transition while (F = l,/i = 0) corresponds to the integrable quantum critical 
point. Throughout this paper, F is set equal to 1 and h is quenched along the dashed line shown with an arrow. 


THE KIBBLE-ZUREK SCALING 

As emphasized in the main text, so far as the slow quenching of the longitudinal field h is concerned (especially, 
for large t), one can equivalently work with the effective integrable Hamiltonian given by 


ileff = 


(S2) 


where 6 is a constant m which is inessential in the argument below, and hence set equal to unity hereafter. Using a 
gauge transformation (which flips the spins of alternate sites) and a duality transformation the Hamiltonian in 
Eq. (S2) can be mapped to an equivalent dual Hamiltonian with a nearest neighbor ferromagnetic (FM) interactions: 
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HeB = 




'i 'i+1 


i 

i i 


(S3) 


which is a FM transverse Ising Hamiltonian in an effective transverse field Fes = \ — h?. We note that H^s with h = 0 
represents a critical Hamiltonian. Using the Fourier transformation followed by the Jordan-Wigner transformation, 
the model can be reduced to a two-level problem in the basis |0) (no fermion state ) and |fc, —k) (a state with a pair 
of fermions with quasi-momenta k and —k, respectively) [ 33111 ]; the reduced 2x2 Hamiltonian is then given by 




—i sin k 

— (1 — h^) -I- cos k 


(S4) 


Analyzing the spectrum , = 2\J {(1 — /i^) — cos k}^ + sin^ fc, it is straightforward to show that the model (S3) 


has three QCPs; the energy gap (2efe) vanishes at critical points at /i = 0 and h = ±-\/2, with the corresponding 
critical wave vector (for which the energy gap vanishes) kc = 0 and tt, respectively. We are however interested in the 


transition at /i = 0 which is the only relevant QCP to the context of the original Hamiltonian (SI). To focus on the 
critical point at /i = 0, we expand the Hamiltonian (p4| in the vicinity of fc = 0 to arrive at the Hamiltonian 


which shows only one quantum critical point at h = 0. Analyzing the simplified form of the spectrum Ck = 
v'(/i2 - fc72)2 one immediately finds for h = 0, the gap {AEk = 2ek) k, yielding z = 1 and for k = 0, 

gap scales as yielding uz = 2, and hence v = 2 (referred to as Vh in the main text). 

Let us now point out that the quenching h = —t/r, with t going from — oo to 0, in the original Hamiltonian is 
equivalent to driving the reduced Hamiltonian ( |S5[ ) from /i —>■ oo to the QCP at = 0 by a non-linear protocol 
(t/r)^; in both the cases the system is initially prepared in its ground state. Even though the non-adiabatic transition 
probability for the mode k (p^) can not be calculated directly using the Landau-Zener formula for such a non-linear 
protocol, one can make appropriate rescaling in the corresponding Schrodinger equations [33113] to argue that it would 
be a function of the dimensional combination of i.e., pk = where E is an unknown scaling function. 

Since the gapless QCP is characterized by gapless excitations k, the scaling of the residual energy can be obtained 
as Cres ~ / dkkJ-{kP‘T'^^^) ^ this matches perfectly with the KZ prediction with d = z = 1 and v = 2 and the 

numerical result presented in the main text. 


SLOW QUENCHING AND NON-ANALYTICITIES IN THE RATE FUNCTION 


We shall now calculate the nature of the Fisher zeros of the effective partition function m obtained from the 
Loschmidt overlap (LO) when the parameter h of the Hamiltonian (S3) is quenched from a large positive to a large 


negative value following the protocol h = —t/r; this is equivalent to the slow quenching of the longitudinal field h in 


the original Hamiltonian (SI). But there is a subtle difference that needs to be emphasized: the field h contributes 
a quadratic h'^ term to the transverse field of the equivalent model (S3), thus as h is linearly changed from a large 


positive value to the negative value in model (SI), the parameter changes from a positive initial value to zero 


(i,e, the QCP) and returns to the original initial value at the final time; this in a sense is a reverse quenching of the 
transverse field of the Hamiltonian (S3) as studied in [ 31 ] in a non-linear fashion. 


To calculate the Loschmidt overlap of a system of length N defined by f{z) = — ln(^/| exp(—iLjz)|'0/)/Af, where 0 
is the complex time, Hj is the final Hamiltonian and ['(/'/) the state reached following the quantum quench, we focus 
on the reduced 2x2 Hamiltonian ( |S5[ ). Summing over the contributions from all the momenta mode, a few lines of 
algebra leads us to the expression [39] 


/(^) = -y ^ln(^(l-pfe)-hpfeexp(-2e{z)) 


(S6) 
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FIG. S2: The rate function for quenching h — 3 to h — —3 shows sharp non-analyticities at periodic intervals in time with 
system size N = 400 and several r’s. 


where pk is the non-adiabatic transition probability for the mode k. The zeros of the “effective” partition function 
(where f{z) is non-analytic) are given by: 


= r7 ( 


Pk 


1 -Pfc 


) + *7r(2n + 1) 


(S7) 


where n = 0,±1,±2,'--. For a non-linear reverse quenching protocol, the expression for pk can not be exactly 
determined using the LZ formula (though an exact form can be obtained for the linear case [9E])- However, it can 
be argued pk = Gi{k — where ko is the wave vector for which pi~ is maximum which shifts to A: = 0 for 

large t and G is an unknown function. We find from Eq. (S7) Fisher zeros cross the imaginary axis for a particular 


value of /c* for which pk^ = 1/2 [321 [10] and the rate function shows sharp non-analyticities at tn = 

For the present case, to calculate the Fisher zeros and especially the rate functions I{t) we shall use the form of the 
Hamiltonian near k = 0 given in Eq. (S5) (to avoid the influence of the QCP at ft. = V^), when ft is quenched from 
-|-3 to —3. Numerically integrating the Schrddinger equation describing the dynamics of the Hamiltonian (S5) with 


the initial condition that the system is in the ground state of the initial Hamiltonian, we obtain the value of pk which 
is then substituted in the expression of the rate function: 


I{t) = “ ^ ^ 4pfc(pfe - 1) sin^ . 


(S8) 


As shown in Eig. |S2[ this qualitatively explains the periodic occurrence of DPTs presented in Fig. [^in the main text. 
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